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Abstract 

Central nuclear collisions at energies far above 1 GeV/nucleon may 
provide for conditions, where the transition from highly excited hadronic 
matter into quark matter or quark-gluon plasma can be probed. Here I re- 
view our current understanding of the physical properties of a quark-gluon 
plasma and review ideas about the nature of, and signals for, the deconfine- 
mcnt transition. I also give a detailed presentation of recent progress in the 
treatment of the formation of a thcrmalized state at the parton level. 

INTRODUCTION 

Overview 

After a collision between a '^^S nucleus of 6.4 TeV total energy and a heavy target 
nucleus several hundred charged particles are emitted. No one in his or her right mind 
would care to study such events, unless there existed a compelling reason for doing 
so. The current interest in nuclear collisions at very high energies (far above 1 GeV/u 
in the cm. system) is fueled by the expectation that a quark gluon plasma may be 
created temporarily in these events^""^. Whereas there is general consensus among 
theorists that QCD at thermodynamic equilibrium exhibits a phase transition from 
the normal color-confined phase of hadronic matter with broken chiral symmetry to a 
deconfined, chirally symmetric phase at sufficiently high energy density, many aspects 
of this transition are still a matter of debate. Such "details" are, e.g., the order of the 
phase transition, the precise value of the critical energy density, the nature of exper- 
imentally observable signatures of the transition, and how fast thermal equilibrium 
conditions are established in nuclear collisions over a sufficiently large space-time vol- 
ume. These questions require much further theoretical and experimental study. Here 
we are concerned with an up-to-date survey of (mostly) theoretical aspects. More 
detailed discussions of selected theoretical topics, as well as reviews of experimental 
results, can be found in other lectures presented at this school. 



The "Cosmic" Connection 

Perhaps the most compelling reason why we should attempt to study the quark- 
gluon plasma transition in laboratory experiments is that it must have occurred in 
the early universe^. The history of the thermal evolution of our universe is depicted 
in Fig. 1. The relation between temperature T and cosmic time t is approximately 
given by^: 

TMeV ^ {5.75Nf{T))h-l^\ 

where Nf{T) describes the number of particle degrees of freedom that act as effectively 

massless modes at a given temperature, and the subscripts indicate the units in which 
T and t are measured. From the present temperature of the cosmic background 
radiation (2.7 K) we extrapolate back to a temperature of about 2 x lO^^K ~ 200 
MeV at about 20/i.s after the "big bang". This is the temperature above which, as 
we will discuss in detail below, hadrons dissolve and their constituent quanta, quarks 
and gluons, are liberated. 
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Figure 1. Thermal history of the Universe. 

Tracing the history of our universe backward, this is only the first phase transition 
involving fundamental quantum fields that we encounter. Most likely, more trans- 
formations of a similar nature have occurred at even earlier times. If our current 
ideas of the origin of electroweak symmetry breaking are correct, a phase transition 
in the Higgs vacmim took place at i w lO^^^s, when the temperature was around 
250 GeV. A similar phase transition at much earlier times {t ~ 10~^^s) associated 
with the "grand" unification of electroweak and strong interactions may have led to 
exponential inflation of cosmic scales, due to the anti-gravitational pressure exerted 
by an unstable vacuum state^'''. 

Why should we care about events occurring at such unimaginably short moments 
after the creation of our universe? The reason is that some of the unsolved problems 
of cosmology are, most likely, associated with events in the era before about 10~^s: 
• The observed baryon number asymmetry in the universe is probably due to baryo- 

gencsis by topologically nontrivial field configurations during the electroweak 

phase transition. 




• The nature of dark matter may be associated with properties of the Higgs or the 
QCD vacuum. 

• The near isotropy of the background radiation, and the deviations from homogene- 
ity imderlying the observed large-scale structure of the universe were, according 
to current thinking, determined during the phase transition causing inflation. 

Understanding of the dynamical nature of phase transitions in elementary quan- 
tum field theories is thus essential to further progress toward the solution of these 
puzzles. Since the QCD phase transition is the only one of these that appears acces- 
sible to laboratory experiments, we must study it with high priority. 

Quark Stars 

In addition to events in the very early universe, quark matter may also play a role 
in the internal structure of collapsed stars. At the high densities reached in the core of 

neutron stars nucleons may well dissolve into their constituents, forming baryon-rich 
cold quark matter. This would not lead to greater stability of neutron stars, quite to 
the contrary: since quark matter would allow for a higher central density of the star 
at a given total star mass, its formation would actually facilitate collapse to a black 
hole. Neutron stars with a quark core have a lower value of their upper mass limit, 
probably somewhere between 1.5 and 2 solar masses®. Stars with a quark core would 
be more compact and hence could sustain higher rotation rates^. This observation 
would be of practical interest, if pulsars with periods in the sub- millisecond range are 
eventually detected. 

The problem with quantitative predictions here is that dense baryon-rich nuclear 
matter is expected to contain a large strangeness fraction, because the inclusion of 
strange quarks can lower the Fermi energy. This holds true for baryonic matter as well 
as quark matter. Unfortunatc^ly the equation of state of baryonic matter containing 
hyperons is poorly known. The scalar coupling strength of A-hyperons, which is 
not well determined experimentally has a large influence on the central density of 
hyperon-rich neutron stars. This uncertainty can easily mask the formation of a quark 
matter core®. Better understanding of the equation of state of hyperon interactions 
in dense baryon-rich matter, which can only come from high-energy nuclear collisions, 
is therefore essential for further progress towards the understanding of neutron star 
structure. 

THE EQUATION OF STATE OF HADRONIC MATTER 

There are two approaches that have been widely followed to the problem of the 
equation of state of strongly interacting matter: 

(a) Consider color-singlet hadrons, i.e. baryons and mesons, only and see how far one 
gets with taking into account their known interactions and excitation spectra. 
This approach was pioneered by R. Hagedorn, and later studied by Hagedorn 
and Rafelski^°, Walecka^^, Gasser and Leutwyler^^, and many others. 

(b) Consider the fundamental constituents of hadrons, i.e. quarks and gluons, and 
study the equation of state on the basis of quantum c:liromodynamics. This 
approach, anticipated in pre-QCD days by P. Carruthers^'^, was pioneered by 
Collins and Perry^^, Baym and Chin^^, McLerran et al.^^ and by Shuryak^''. 

Here we will start by examining the approach based on hadrons and their in- 
teractions, because its results are important in their own right. Any quark-gluon 
plasma formed in a nuclear reaction will eventually hadronize and evolve through a 



hadron-dominatcd break-up phase. In addition, potcintial signaturc;s for qnark-ginon 
plasma formation must always be compared with predictions of a scenario that makes 
no reference to a phase transition from color-singlet hadrons to deconfined quarks and 
gluons. 

Hot and Dense Hadronic Matter 

How far do wc get in explaining the equation of state of hadronic matter by con- 
sidering hadrons alone? Experimental observations of baryon and meson resonances 
indicate that the mass density of hadronic states grows exponentially, 



as originally postulated by Hagcdorn^*^. It is now understood that an exponential 
mass spectrum is the natural consequence of quark confinement, and it is found in 
even the simplest hadron models that incorporate the confinement concept, such as 
the string model^'' or the MIT-bag model^*^. A numerical simulation of the spectrum 
of excited string modes agrees very nicely with the observed hadron spectrum even 
at low energies^ ^. 

Now consider highly excited hadronic matter, characterized by a temperature T. 
The energy density of excited states is then given by (2.1) integrated over momentum 
space, i.e. 



Obviously, n{E) is only integrable as long as the factor in the exponent remains 
positive, i.e. for T < Tc = mo. From experimental data we know that toq ~ 
200 MeV, i.e. the temperature of a (noninteracting) hadronic resonance gas cannot 
exceed the Hagedorn temperature Tc ^ 200 MeV. Cabibbo and Parisi pointed out 
in 1975 that the distribution (2) remains integrable at T = Tc when a < —4, and 
then corresponds to a finite energy density £ = / EdE n{E) at the critical point^^. 
This is reminiscent of the singular behavior of thermodynamical quantities at a first- 
order phase transition, and they speculated that it signals the transition to a color- 
deconfined quark-gluon plasma. 

Although Hagedorn has stressed that the inclusion of hadronic resonances eff'ec- 
tively accounts for a large part of the interactions among hadrons, not all interactions 
are taken into account in this simple manner, especially if one neglects the resonance 
widths. When Hagedorn and Rafclski^'' studied the effect of the finite internal size of 
hadrons by utilizing the finite excluded volume approximation, they found that the 
equation of state developed a singularity at finite energy density at a slightly lower 
critical temperature Tc < mo, independent of other details of the hadronic mass spec- 
trum (2.1). We conclude from all this that hadronic matter as composed of individual 
color-singlet hadrons ceases to exist at temperatures exceeding Tc ~ 150-200 MeV. 
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Medium Modifications of Hadron Structure 



It is useful to sharpen the question concerning the influence of interactions on 
the hacironic equation of state and ask whether hadron masses themselves depend 
on temperature and density. As this interesting problem has been approached from 
several different angles, it is quite enlightening to discuss a few of these. 

(a) Sigma- condensate models: 

In these models the nucleon mass M is generated by the coupling to a scalar field a: 

jCn = ipi-y ■ dtp — gtptpa, (3) 
which develops a nonvanishing vacuum expectation value 

ao = (0|a|0) = M/g. (4) 
However, since the cr- field interacts with hadrons, via an equation of the form 

id''+til){a-ao) = -gm, (5) 
its expectation value is shifted in dense nuclear matter according to the equation 



Here 



M* =ga = M- 
d^k M* 



{n{k)+n{k)) 



(6) 



(7) 



(27r)3 E{k) 

is the scalar baryon density, E{k)'^ = fc^ + M and n{k),ri[k) denote the momentum 
space density of baryons and antibaryons, respectively. The numerical evaluation of 
eq. (2.6) shows that M* drops rather suddenly to a very small value at T ~ 200 
MeV at zero net baryon density p^^ . The decrease with p at T = is more gradual, 
already corresponding to M* /M ~ 0.6 at nuclear matter saturation density po- 
becomes very small at p/po ^ 3 — 4. 





Figure 2. EfTcctivc mass of the nucleon in the a—uj model: (a) as function temperature T at p=Q, 
(b) as function of the baryon density p at T = 0. 



(b) Dispersion relations: 



A second method of approach makes use of available data on hadron-hadron 
scattering via dispersion relations. The propagator of a hadron, say a pion or kaon, 
in medium is modified due to interactions with other particles. In many cases this 
interaction is dominated by resonance scattering, schematically illustrated in Figure 
3. This is the case, e.g., for pions which interact with other pions from the medium 
to form a short-lived rho-meson, and with kaons that form a K* resonance. This 
also occurs with mesons interacting with nucleons via the S* resonance at 1385 
MeV, but not with mesons which have a dominantly nonresonant interaction with 
nucleons. 
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Figure 3. Resonance scattering on pions from the medium leads to a medium-dependent mass shift 
of (a) pions, (b) kaons. 

The in-medium propagator of a (pseudo-) scalar meson has the form 

D{k) = - - ni(fc))-i = {k^ - m*(fc)2)-i (8) 

where i denotes the meson species and the polarization function Tli{k) can be ex- 
pressed as the sum of the contributions from interactions with all other hadrons 
contained in the medium^^: 

Here fij{0, s) is the forward scattering amplitude at energy s = {p + k)"^ and nj{p) 
denotes the momentum space density of hadrons of species j in the medium. Numer- 
ical evaluation of expression (8) for pions and kaons reveals rather small mass shifts 
Ami = = 0) — TOj of order 10-30 MeV in baryon-symmetric hadronic matter 

up to temperatures around 150 MeV, beyond which the approximations entering into 
(9) become unreliable. 

However, it turns out that the effect on kaons in baryon-rich matter can be 
substantial. As shown in Figure 4, the mass shift of kaons at T = 160 MeV and 
baryon density p = Apo is of the order of 100 MeV, being positive for and negative 
for K~ . The reason for the opposite behavior is that the interaction of JsT" -mesons 
with nucleons is resonance-dominated and attractive, whereas the interaction of K^- 
mesons is nonresonant and repulsive. An immediate consequence of the different mass 
shift is that the relative abundance of charged kaons in dense baryon-rich matter is 
modified^'''. In fact, the mass shift acts counter to the effect of the nonvanishing 
baryochemical potential and limits the ratio of abundances K~^/K~ to about 5 at 
T c± 150 MeV even for extreme baryon density. Notably this is just the ratio observed 
in heavy ion collisions at the AGS^^. Is this a coincidence? Although kaons in these 
reactions are probably emitted long after the moment of highest density in the fireball, 
the K~^/K~ ratio may be determined by the reactions occurring at or before the 
moment of highest compression. The abundances probably get out of equilibrium 
and remain frozen as the fireball expands and breaks up rapidly. 
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Figure 4. Mass shift of kaons in baryon-rich hadronic matter. The figure shows the medium- 
dependent change in the energy of positive kaons (soUd lines) and negative kaons (dashed hnes) as 
function of momentum for T = 160 MeV and twodifferent baryon densities. 

(c) QCD sum rules: 

A third avenue of approach to medium modifications of hadronic properties 
touches base with the underlying QCD dynamics by identifying the sigma conden- 
sate (To with the scalar quark- antiquark condensate {qq) in QCD. This permits one to 
advance the argument^^ that all hadron masses should be modified in the same way, 
according to 

m* (med|gg|med) 



rrii 



lw|o> 



(10) 



A qualitative estimate of the change in hadron masses can then be obtained from any 
estimate of the scalar quark density in the medium, as compared with the vacuum 

condensate^® 

{0\qq\0) = -(225 ± 25MeV)^. (11) 

A simple way to estimate the medium contribution to {qq) is to calculate the scalar 
quark density contained in the thermal pions 



5{qq)-K 



(2^ 



{TT{k)\qq\TT{k))n^{k), 



or, for baryon-rich matter, in the nucleons 

d^k 



5{qq)N j 



{N{k)\qq\N{k))nN{k). 



(12a) 



(126) 



Eq. (12a) may be conveniently evaluated in the framework of chiral perturbation 
theory Results are shown in Figure 5a. More model- independent estimates for 
(12a,b) can be obtained by the so-called QCD sum-rule approach^^"^^. Figure 5b 
shows the predictions of this technique for the medium modifications of the conden- 
sates of light quarks, strange quarks and gluons in the QCD vacuum as function of 



baryon density (at T = 0) . Obviously, the condensate of light quarks is most strongly 
affected in both approaches, indicating a dramatic change in hadron structure at 
200 MeV or /9 ~ 4po- 




T p/Po 

Figure 5. Medium modification of QCD condensates, (a) Change of the quark condensate in a pion 
gas as function of temperature^^ . (b) Dependence of the condensates of Ught quarks, gluons and 
strange quarks as functions of the baryon density"^^ . 

Overall, the agreement between the predictions of the different approaches agree 
remarkably well, reinforcing faith in their accuracy. In general, one finds that hadron 
structure is much less disturbed by a baryon symmetric medium of finite temperature 
than by a medium with finite net baryon density. The right place to look for medium 
effects on hadron structure is therefore in relativistic nuclear collisions involving a 
large degree of baryon stopping, such as the Brookhaven AGS especially with the 
^^''Au beam — or possibly the CERN-SPS in Pb+Pb collisions. In the case of matter 
depleted of baryons the medium effects set in so late that it will be very difficult to 
disentangle them from the quark-gluon plasma phase transition. 



QUARK-GLUON PLASMA 

The theory of the equation of state of quark matter is conceptually much simpler, 
because it is directly based on the fundamental QCD Lagrangian 

Nf 

CQCB = -\J2K^P''"' + Jl^('^'^^-9rA;^-mf)i; (13) 

a /=! 

where the subscript / denotes the various quark flavors u,d,s,c, etc., and the non- 
linear glue field strength is given by 

f;, = d^Al - d^Al + gfabcAlAl. (14) 

QCD predicts a weakening of the quark-quark interaction at short distances (or high 
momenta Q^), because the one-loop series for the gluon propagator yields a running 
coupling constant 



The QCD stale; parameter is now quite well determinecP^ to be A ~ 150 MeV. 
The vanishing of the QCD coupling constant at short distances, called "asymptotic 
freedom" , has often been taken to imply that interactions among quarks and gluons 
arc negligible in the limit of high temperature or high density. As we shall sec below 
this is not the whole truth, because long-wavelength modes continue to have an 
important influence on the properties of the quark-gluon plasma. 



Equation of State of Quark Matter 

So let us first suppose that interactions among quarks and gluons are negligible 
at high energy density and see what we get. At temperature T and quark chemical 
potential fi (one-third of the baryochemical potential /xb = 3/i), the energy density of 
free gluons, quarks and antiquarks is (a detailed derivation can be found in ref. 1): 

Eg = 16^T^ (16a) 
e,+e^= 6iV, (I^T^ + \^?T- + . (166) 

Since we have neglected the quark mass, we inserted a factor Nf counting the number 

of quark flavors active at a given temperature (essentially those with mj < T). 
Clearly Nf > 2 at all relevant temperatures, and Nf = 3 in the range 200 MeV 
<^ T < 1 GeV where strange quarks can be considered as light, as well. The other 
thermodynamical quantities of interest, i.e. pressure P, entropy density s, and baryon 
number density pB, are easily obtained from eqs. (16): 

Note that the simple relation between e and P holds only as long as all particles are 
considered massless. The general relation is 

In order to find the location of the phase transition toward hadronic matter we have 
to incorporate the breaking of scale invariance provided by QCD interactions. The 
simplest way of achieving this is by invoking the argument associated with the MIT- 
bag model: Free quarks and gluons can only propagate where the complex structure of 
the real QCD vacuum has been destroyed. The vacuum realignment costs a certain 
amount of energy per unit volume, expressed by the MIT-bag constant Eq = B 
(150-200 MeV)^, and the "wrong" vacuum is endowed with a negative pressure 
Pq = —B. The relation eq = —Pq is mandated by the Lorentz invariance of the 
vacuum state. The negative sign is easily understood as signal of the instability of 
the wrong vacuum state which collapsc;s if not supported by the pressure; provided 
by free partons propagating in the volume filled with it. The equation of state of the 
quark-gluon plasma then takes the simple form 

E = Eg{T, H) + £g{T, /i) + £g{T, fl) + B' (Wtt) 

P = Pg{T, n) + Pg{T, p) + P^(r, p) - B. (196) 



A lower limit of stability of the plasma state is obtained by setting P = 0. This yields 
the stability line in the T-fi plane shown in Figure 6. Of course, the plasma phase 
becomes unstable against formation of a gas of color-singlet hadrons even earlier, 
when its pressure equals that of a hadron gas at the same temperature T and chemical 
potential fx. 
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Figure 6. Stability line of the quark-gluon plasma in the T-/i plane for two values of the strong 
coupling constant. The lines indicate where the pressure of a gas of quarks and gluons vanishes. 
The quark-gluon plasma is unstable in the lower left region of the figure. 

In order to explore this aspect further, let us restrict ourselves to the baryon-frce 
case /X = 0, and explore the coexistence between the quark-gluon plasma phase given 
by eq. (19a) and hadronic matter, represented by a gas of massless, noninteracting 
pions: 

e^ = 3^T\ P^ = y^. (20) 
e, P, e„ and Pj^ are shown as functions of T in Figure 7. 
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Figure 7. Energy density and pressure of a pion gas and a free quark-gluon gas as function of 
temperature. A first-order phase transition occurs where the pressure of the two phases is equal. 
The solid line indicates the stable phase. 



Thcrniodynamical phase; stability requires that the phase with the larger pressure 
dominates, and phase equilibrium is achieved when -P(T'c) = Pti{Tc)- As Figure 7b 
shows, one finds ~ 0.75Bi ~ 150 MeV with B ~ (200 MeV)'') in this model. 
Due to the vacuum rearrangement energy B, the energy density between the two 
phases differs greatly at this point, by the amount A^c — ~ 0.8 GeV/fm'^. This 
simple model obviously predicts a first-order phase transition between the pion gas 
and quark-gluon plasma with a large latent heat Aec. 

Of course, our model is grossly oversimplified because, as we saw earlier, other 
hadron masses begin to decrease substantially around T ~ 150 MeV, leading to an 
increase in the energy density of the hadronic phase. In parallel, interaction between 
quarks and gluons cause a reduction in the energy density and pressure of the plasma 
phase. In first-order perturbation theory, the modification of the plasma equation of 
state is (for Nf = 2): 
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yielding a reduction by about a factor 2 when = 0.5. 

More reliable predictions concerning this phase transition can at present only 
be obtained by numerical simulations of the QCD equation of state on a discretized 
volume of space-time, usually referred to as lattice gauge theory. In this approach^^ 
one approximately calculates the partition function for a discretized version of the 
QCD Lagrangian (13) by Monte-Carlo methods. In principle, this technique should 
accurately describe the quark-gluon plasma as well as the hadronic phase but, in 
practice, its accuracy especially at low temperature is severely limited by finite size 
effects and other technical difficulties. Where the numerical results are most reliable, 
i.e. for the pure gluon theory without dynamical quarks, the calculations predict a 
sudden jump in the energy density at a certain temperature while the pressure rises 
more gradually, as shown in Figure 8. 



CQ 

m 



1 

.8 
0.6 
0.4 
0.2 




J I \ L_ 



_J \ I L_ 



□ £ 



SU(3) Lattice 



1.5 2 



2.5 



Figure 8. Energy density e and pressure P of pure glue matter as calculated by simulations of SU(3) 
lattice gauge theory^'', e and P are plotted relative to the Stefan-Boltzmann limit, eq. (16a). 

When dynamical quarks are added, the picture becomes less clear for two reasons. 
One is that the calculations involving fermion fields on the lattice are much more 



time consuming, and hence tlie numerical results are less statistically meaningful and 
reliable. Moreover, the definition of quark confinement becomes rather fuzzy in the 
presence of light quarks, because the color flux tube between two heavy quarks can 
break by creation of a light quark pair: QQ {Qq){(lQ)- E-g-, highly excited states 
of charmonium can break up into a pair of D-mesons. Thus, in the calculations the 
QQ potential does not rise linearly with distance, but is effectively screened. 

For masslcss dynamical quarks there exists a new order parameter, the quark- 
antiquark condensate (0|gg|0) in the vacuum. When it assumes a nonzero value, 
chiral symmetry is spontaneously broken, as can be seen as follows: The scalar quark 
density has the chiral decomposition qq = qLQR'^^^n.QL j hence the broken vacimm state 
contains pairs of quarks of opposite chirality. A left-handed quark, say, can therefore 
annihilate on a left-handed antiquark in the vacuum condensate, liberating its right- 
handed partner. This process is perceived as change of chirality of a free quark, which 
is exactly the same result as that of a nonvanishing quark mass. However, in reality 
the mass of the light quarks u, d is nonzero, and the chirality of a light quark is never 
exactly conserved. 

All one can do, therefore, is to look for sudden changes in the distance at which 
color forces are screened, or in the quark condensate. If these are discontinuous, 
one deals with a phase transition, otherwise with a possibly rapid, but continuous 
change of internal structure as it occurs, e.g., in the transformation of an atomic gas 
into an electromagnetic plasma. The identification of the nature of the phase change 
is complicated by finite size effects. The best published results, by the Columbia 
group for a 16"^ X 4 lattice indicate a surprisingly strong dependence of the phase 
diagram on the magnitude of the strange quark mass. For the physical mass rUs — 
150 MeV there seems to be no discontinuity, but only a rapid change in the energy 
density over a small temperature range (about 10 MeV). However, it is probably 
premature to consider this as the final word. 

Intermezzo: Astrophysical Implications 

We will return to the results of lattice gauge theory in a moment, but let us 
pause briefly to study some astrophysical implications of a hadron-quark-gluon phase 
transition. If it is indeed of first order, as the results discussed so far suggest, both 
phases can coexist over a certain range of temperatures. In particular, the quark- 
gluon phase can exist in a supercooled state between the critical temperature Tc, 
where its pressure falls below that of an isothermal hadron gas, and the temperature 
T, where its pressure becomes negative signaling intrinsic instability of this phase 
(see Figure 7b). The conversion of the metastable quark-gluon phase into hadrons 
then proceeds by the formation and subsequent growth of bubbles of hadronic gas 
imbedded in the plasma. The dynamics of bubble formation is governed by the free 
energy difference as a function of bubble radius R: 

AE{R) = -AP—R^ + a4TrR^+aR+..., (22a) 
o 

where AP > is the difference in pressure between the two phases, a is the surface 
tension of the interface, and a is a coefficient related to the Casimir energy of the 
bubble configuration. The function AE{R) is schematically depicted in Figure 9a. 
The notable feature here is the fact that bubble formation is thermodynamically dis- 
couraged for small bubbles. Only when a bubble grows to a radius R > Rc, where Rc 
denotes the location of the maximum AEc 2a/ AP, due to some thermal fluctua- 
tion will it continue to grow on its own. As the probability of forming such a critical 
droplefi^ is obviously hindered by a factor exp(— AiSg/T), the suppression depends 



scnsitivcily on the size of the interface energy a. Its value has recently been calculated 
in the framework of QCD-related models^^ as well as by lattice simulations"^''. All 
results indicate that the surface tension should be small: 



a ~ 20-50 MeV/fm^ 



(23) 



implying that it is probably difficult to supercool a quark-gluon plasma for a consid- 
erable time (see ref. 38 for a recent discussion). 





Figure 9. Energy of a hadronic bubble in tlie quark-gluon plasma versus its radius: (a) for T<Tc, 
(b) for T>To. Rc in (a) indicates the critical droplet. The solid curve in (b) predicts the "Swiss 
cheese" instability of a quark-gluon plasma against formation of small inhomogeneities. 

On the other hand the small surface tension raises the importance of the next 
term in the series (22a), which is linear in the bubble radius R. That may have a 
dramatic effect when we consider the opposite case, namely, superheated hadronic 
matter droplets forming inside a large quark-gluon plasma bubble. This is precisely 
the situation one hopes to find right after a highly energetic nuclear collision. When 
we now ask about the stability of the plasma, we must consider the reverse situation 
described by a droplet energy obtained from (22a) by altering the signs of all terms 
with an odd power of R (see solid curve in Figure 9b): 

4:7T 

AE'{R) = AP—R^ + gA-kR^ -aR+.... (226) 
o 

The change in sign of the linear (as well as the cubic) term was pointed out by 
Lana and Svetitsky"^^ who went on to study its possible consequences. Since the 
linear term dominates AE'{R) < for small R, and it is favorable to break up the 
plasma by forming small droplets of superheated hadronic matter inside. This has 
been termed the "Swiss cheese instability" of the quark-gluon plasma. More detailed 
studies show the hadronic droplets to be unstable against deformations, pointing to 
the possibility that the plasma phase may be characterized by a complex structure 
of inhomogeneities in the region of phase coexistence. There is preliminary support 
for this exotic picture from lattice simulations^*^. 

Perhaps the most important effect caused by coexistence of the; two phases is 
due to the difference between the baryon densities in them at equilibrium'*^. As 
discussed above, the baryon density in the quark-gluon phase at given temperature 
T and baryochemical potential fiB is (using eqs. 21, 17): 

Pi^' - 1^ {t^ + ^^l) , (24) 



whereas the net baryon density in the hadronic phase (counting only nucleons, not 
excited baryon resonances): 

P^f^-8?(^)'''e-'"~/^ (25) 

where mjv is the nucleon mass. The ratio / p^'' is very sensitive to the precise 

value of the critical temperature T^, but generally one finds p'^^ » , as illustrated 
in Figure 10. 

Consequently, if the transition proceeds at or near chemical equilibrium, baryons 
become enriched in the quark-gluon plasma phase. For some time it was speculated 
that this effect could result in local fluctuations of the proton-neutron ratio in the 
early universe that were large enough to influence the course of nucleosynthesis. The 
argument goes roughly as follows: Neutrons will quickly diffuse away from a local 
surge in the baryon density, leaving behind a region of abnormally large proton- 
neutron ratio. Nucleosynthesis there will result in a reduced number of neutron-rich 
isotopes, such as ^He and ^Li. However, if the size of the region of anomalous baryon 
concentration is too small (less than about 20 m) neutrons can diffuse back during the 
era of nucleosynthesis and even things out. For the small value of the interface tension 
(T, eq. (23), the effect of local inhomogeneities in the baryon density is probably 
negligible^^'^^. 




Mb (MeV) 

Figure 10. Ratio of the baryon densities of the quark-gluon plasma and a hadronic gas at the 
coexistence line as function of the bajryochemical potential /hb- Conditions in the early universe 
correspond to the line /iB=0. 

In rclativistic heavy ion collisions, the most important result of the concentration 
of baryons in the quark matter phase may be that it is accompanied by an enrichment 
with strange quarks (as opposed to strange antiquarks). The interest in such a mech- 
anism derives from speculations that quark matter droplets with a high strangeness 
content (so-called strangelets) might be stable or metastable**^'**^, because a large frac- 
tion of strange quarks allows to construct electrically almost neutral nuclei with large 
baryon number without raising the proton- neutron asymmetry. As Liu and Shaw'*^, 
as well as Greiner, Koch and Stocker^'' showed, the quark matter — hadron gas phase 



transition effectively acts as a distillation process for strange quarks, if there is a 
net baryon surplus'**. Several experiments at Brookhaven (AGS) and CERN (SPS) 
are searching, or preparing to search, for strangelets produced in nuclear collisions*®. 
Their existence would have far-reaching technological implications^'^. 

INTERACTING QUARK-GLUON PLASMA 

Let us return now to our main theme, i.e. the physical properties of a quark- 
gluon plasma. So far, we have neglected interactions among quarks and gluons in the 
deconfined phase, except in eq. (21), where we included the contributions of one- loop 
diagrams to the energy density and pressure of a quark gluon gas: 

V 47r V 15 V 2l7r J 10 

+ {'-^2-s)^,^{n^T^ + l,^), (26) 

where we have omitted the vacuum energy B, because it is a non-perturbative con- 
tribution. Graphically, eq. (26) corresponds to the diagrams 





(27) 



where dashed lines denote gluons, and straight lines denote quarks. As (26) shows, we 
still have P = |e. This changes in the next order, because the two-loop contribution 
to the gluon energy density is found to diverge! The physical reason for this divergence 
is that gluon and quark degrees of freedom develop an effective mass, which leads 
to screening of long-range color-electric forces. Technically, the screening mass is 
obtained by summing an infinite chain of one-loop insertions in the gluon propagator 

+ --\^) + +■■■ (28) 



The contribution of all diagrams except the first two (which are already included 
in (26)) can be summed analytically and yields a contribution to the gluon energy of 
order a"^/-^ with a rather large coefficient^^ . 

The QCD Plasmon 

We can obtain more insight into the properties of the interacting quark-gluon 
plasma by looking at the gluon propagator, represented graphically in (28), itself. 
Because of gauge invariance, k^Dn^ik) = 0, it can be decomposed into a longitudinal 
and a transverse part, which are scalar functions of the variables w = k^ and k = |k|. 
These arc most conveniently written in a form borrowed from electrodynamics of 
continuous media: 

^-(-'^) = Mc^,fe)l^-fe^ ' ^'''^ 



where the color-dielectric functions are given by^^: 

(jj + k 



£z,(a;, fc) = 1 
ST{co,k) = 1 
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Several things arc noteworthy about eqs. (29, 30). First they imply that static 
longitudinal color fields are screened: 



Dl{0, k) 



1 



ei(0,fc)A:2 k^+g^T^ fc2 + A^' 



2 • 



(31a) 



The Debye length obviously is Ad = {gT)~^ . On the other hand, cqs. (29b, 30b) 
show that static transverse (magnetic) color fields remain unscreened at this level of 
approximation: 

£>T(0,fc)=-^. (316) 

The static magnetic screeening length is of higher order in the coupling constant; 
lattice gauge calculations^^ have shown that X'^ = Cg'^T. 




Figure 11. Gluon dispersion relation in the perturbative quark-gluon plasma. All quantities are 
measured in units of the effective gluon mass m*=gT/\/3. Solid line: longitudinal plasmon mode; 
dashed line: transverse collective gluon mode. 

For a finite frequency w the in-medium propagators (29) have poles corresponding 
to propagating, collective modes of the glue field. The dispersion relation for the 
longitudinal mode: 

ez,(a;,fc) = 0, (32a) 



called the plasmon, has no counterpart outside the medium. The analogous relation 
for the transverse mode: 

eT{co,k) = k^/u^ (326) 



describees the; ciffccts of the medium on the free ghion. The bcihavior of both modes is 
remarkably similar. For k ^ they yield an effective gluon/plasmon mass 

Wi,WT^ m; = -^ffT, (33) 

whereas for large momenta {k ^ oo) one finds 



UL{k)^k, a;T(fc)^ ^fc2 + -<?2T2. (34) 

The full dispersion relations are shown in Figure i 1 . 

For plasma conditions realistically attainable in nuclear collisions (T ~ 250 MeV, 
g ~ ^/Anag ~ 2) the effective gluon mass is m* ~ 300 MeV. We must conclude, 
therefore, that the notion of almost free gluons (and quarks) in the high-temperature 
phase of QCD is quite far from the truth. Certainly, one has m* ^ T when g ^ 1, 
but this condition is never really satisfied in QCD, because g — \ even at the Planck 
scale (lO^^GeV), and g < 1 only at energies above 100 GeV. Let us discuss some 
consequences of these results: 




(fm) 



Figure 12. Effective quark-antiquark potential in QCD. Solid line: confining potential of a free QQ 
pair; dashed line; screened potential, eq. (35), of a QQ pair imbedded in the quark-gluon plasma. 

(1) The potential between two static color charges, such as two heavy quarks, is 
screened in the quark-gluon plasma phase. The Fourier transform of eq. (31a) 
yields the potential 

^qq(^) ^ y^"^'' (35) 

with screening length \d ~ 0.4 fm at T = 250 MeV. This screened potential 

is compared in Figure 12 with the confining potential between a free quark- 
antiquark pair. This screening of long-range color forces is, of course, the origin of 
quark deconfinement in the high-temperature phase. An important consequence, 
to be discussed later in the section on plasma signatures, is the disappearance of 
the bound states of a charmed quark pair (cc) in the quark-gluon plasma^^. 



(2) The color screening at large distances cures most infrared divergences in scatter- 
ing processes between quarks and gluons. A self-consistent scheme implementing 
this mechanism has been devised by Braaten and Pisarski^^. It involves the re- 
summation of gluon loops involving gluons with momenta of order gT and has 
been shown to be gauge invariant when also vertex corrections are taken into 
account. 

(3) The finite effective gluon mass m* leads to the suppression of long-wavelength 
gluon modes with k < gT in the quark- gluon plasma. As a result, the relation 
P = is violated and the pressure is reduced. We now have two mechanisms 
that can be responsible for P < ^e: an effective gluon mass m* and a nonva- 
nishing vacuum energy B. A fit to recent SU(3) lattice gauge theory results^^ 
with m* and B taken as free parameters shows that probably both mechanisms 
are at work^'' (see Figure 13). 
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Figure 13. Effective mass m* and vacuum energy density B as function of temperature (lower part) 
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as deduced from results of latt 
P of a gluon gas (see Figure 8). 



57 56 
as deduced from results of lattice gauge theory simulations for the energy density e and pressure 



Properties of the Quark-Gluon Plasma 



(a) Thermalization: 

For those of us interested in the detection of a quark-gluon plasma in nuclear 
collisions it is imperative to know something about its rate of thermalization. Does 
it thermalize sufficiently fast, so that a thermodynamical description makes sense? 
In the picture based on quasi-free quarks and gluons moving through the plasma, 
thermalization proceeds mainly via two-body collisions, where the color force between 
the colliding particles is screened, as illustrated in Figure 14. The technically easiest 
way of looking at this is to consider the quark (gluon) damping rate 7. i.e. the 
imaginary part of the quark (gluon) self energy. For quarks and gluons of typical 
thermal momenta {p ~ T) one finds with the help fo the techniques discussed in the 
previous section^^ 



7, = -4^Im[Tr(7.pI](p))] 
7g = -lm[uT{p)] ^ llg, 



\a,T 1 -Fin ■ 



(36a) 
(366) 



where the factor | reflects the ratio of the quadratic Casiinir invariants C2 of the 
octet and triplet representations of color-SU(3). Gluons simply have a higher color 
charge than quarks (C2 = 3 versus C2 = |) and therefore scatter more often. 

1 1 1 

4 4 4 
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Figure 14. QCD diagrams describing scattering processes contributing to thermalization of the 
quark-gluon plasma. 

One can argue that a better way to look at thermalization is to consider the rate 
of momentum transfer between particles, i.e. weighting the scattering diagrams of 
Figure 14 by sin^ 9, where 9 is the scattering angle. Here one finds^^: 

rf'^ ~ 2.3alT In ^ ; T^*""^ ~ Srf ^ . (37) 

The transport rate F^*''^ is closely related to the shear viscosity of the quark gluon 
plasma. Both approaches yield quite similar numbers for values of the strong coupling 
constant in the range = 0.2-0.5. The characteristic equilibration time, as defined 
as the inverse of the rates 7j or Fj are: 

Tg ~ 1 fm/c, T, ~ 3 fm/c, (38) 

i.e. gluons thermalize about several times faster than quarks^^. Initially, therefore, 
probably a rather pure glue plasma is formed in heavy ion collisions, which then 
gradually evolves into a chemically equilibrated quark-gluon plasma. We shall later 
return to this point in the context of our discussion of microscopic models of relativis- 
tic nuclear collisions, but it is worthwhile to point out one consequence: During its 
hottest phase the QCD plasma is mainly composed of gluons, which are not accessible 
to electromagnetic probes, such as lepton pairs and photons. The gluon plasma can 
only be probed by strongly interacting signals, such as charmed quarks or jets. 

(b) Energy loss of a fast parton: 

One possible way of probing the color structure of QCD matter is by the en- 
ergy loss of a fast parton (quark or gluon). The mechanisms are similar to those 
responsible for the electromagnetic energy loss of a fast charged particle in matter, 
i.e. energy may be lost cither by excitation of the penetrated medium or by radiation. 
Although radiation is a very efficient energy loss mechanism for relativistic particles, 
it is strongly suppressed in a dense medium by the Landau- Pomeranchuk effect^°. 
In the case of QCD this effect has recently been analyzed comprehensively^-'^, and 
the suppression of soft radiation is now firmly established. This limits the radiative 
energy loss to about 1 GeV/fm. Excitational energy loss proceeds via collisions with 
quarks and gluons from the plasma. Here again color screening plays an essential 
role, reducing the rate of energy loss to about 0.3 GeV/fm for a fast quark^^'^^. 



Entropy production 



The correspondence principle asserts that highly excited states of a quantum sys- 
tem usually exhibit quasi-classical behavior. Similarly, the high-temperature limit of 
a quantum system is quasi-classical, since thermal fluctuations dominate over quan- 
tum fluctuations.* Since the quark-gluon plasms is the high-temperature phase of 
QCD, one may ask whether classical, thermal QCD can yield useful results. In view 
of the remark made in the footnote, most predictions of classical thermal QCD will 
be cut-off dependent and therefore not completely reliable, except possibly for those 
quantities that do not involve Planck's constant h. A simple dimensional analysis 
reveals that the combination g'^T does not contain factors of U, and therefore one 
may suspect that the thcrmalization rate 7^, eq. (36b), can be obtained by classical 
considerations. This permits an entirely different, nonperturbative determination of 
the thermalization rate of gluonic matter by simulation of the dynamics of SU(3) 
gauge theory in real time on a lattice. 

The approach to thermal equilibrium in a classical dynamical system is governed 
by the rate of entropy production. More precisely, thermodynamic (ensemble) aver- 
ages can be applied to an isolated classical system, if they coincide with the long-time 
averages, e.g. of a quantity A{t): 

^J^dt A{t) (A)ensemble; (39) 

such systems are called ergodic. It was shown in the important work of Krylov, 
Kolmogorov, and others^^ that ergodicity is really of practical use only when the 
system exhibits the property of mixing, meaning that the equality (39) is approached 
uniformly throughout almost the entire phase space at an exponential rate. This 
condition, in turn, is satisfied if classical trajectories x{t) are unstable against small 
fluctuations almost everywhere in phase space. Such systems are also called strongly 
chaotic. 




Figure 15. (a) Stable phase space trajectory; (b) unstable trajectory. 

Consider two neighboring trajectories, i.e. solutions of the classical equations of 
motion, x(t) and x'(^) = x(i) -|- 6x{t). The trajectory is unstable, if the norm of an 
infinitesimal derivation grows exponentially with time 

\5x{t)\ Doexp{Xt), A>0, (40) 

* As is well-known from the case of electromagnetic black-body radiation, the 
classical limit is only reached for long-wavelength modes, but not for modes with 
wavelength A <C h/T. The quasi-classical description therefore requires an ultraviolet 
cut-off, eliminating these modes. 



as illustrated in Figure 15. A is called a Lyapunov exponent. In general, chaotic sys- 
tems possess several positive Lyapunov exponents, depending on the direction of the 
initial fluctuation (5x(0), but usually the maximal Lyapunov exponent Aq dominates 
in eq. (40) for any arbitrarily chosen (5x(0). 

We are now ready to discuss entropy growth. For a classical system, entropy is 
defined in terms of the volume in phase space covered by an ensemble of identical 
systems:* 

S = In(Ar) + const. (41) 

At first one might think that AF remains constant in time for a Hamiltonian system 
on account of Liouville's theorem. However, for a strongly chaotic system this notion 
conflicts with the finite resolution of any measurement, ultimately with the quantum 
uncertainty limit . The occupied phase space volume AF must, therefore, be 
smeared out with the finite resolution: 

AF^AF, S' = ln(Ar)+const., (42) 

a process usually referred to as coarse- graining. For a strongly chaotic, or mixing 
system the instability of trajectories leads to a growing filamentation of the phase 
space volume occupied by an ensemble, as illustrated in Figure 16. 



Figure 16. Filamentation of the oecupied phase space volume permits a chaotie system to "fill" the 
available phase space homogeneously without violating Liouville's theorem (from Zaslavsky^^). 

A careful analysis^^ of this process shows that the coarse-grained entropy of a strongly 
chaotic system grows linearly with time: 

ds/dt = {h}r, (43) 

where the right-hand side denotes the phase-space average of the sum of all positive 
Lyapunov exponents: 

/i = ^A„^(A„). (44) 

In this way entropy growth, and ultimately thermalization, in classical dynamical 
systems is intimately connected to the instability of its trajectories in phase space. 

One may ask whether this has any relevance for quantum systems? Intuitively, 
one would suppose that it must, because highly excited states of a quantum system 

* The constant is undefined in classical physics; its value is fixed when the volume 
AF is measured in units of , where N is the number of degrees of freedom. 



usually (^xliibit many aspects of classical (lyiiainics. However, the notion of a Lya- 
punov exponent has no immediately counterpart in quantum mechanics, because one 
cannot define a trajectory in phase space for quantum mechanics. Recently, how- 
ever. Gun He has pointed out that there exists a closely related concept in quantum 
mechanics, namely the growth of local inhomogeneities in the Wigner function of an 
ensemble of quantum systems^^. In simple, analytically tractable examples one finds 
that the exponential growth of the variance of this generalized phase space distribu- 
tion is governed by the same exponent as the divergence of trajectories in the classical 
system. This can be utilized to introduce a meaningful concept of entropy for quan- 
tum systems that differs from von Neumann's information entropy, but provides a 
useful description of the approach to equilibrium^^. 

Now back to QCD! It has been known for a long time that non-abelian gauge 
theories exhibit elements of chaos in the classical limit^^. Recently, Trayanov and 
myself showed^^ that randomly chosen field configurations in SU(2) lattice gauge 
theory are characterized by a universal Lyapunov exponent, which scales with the 
average energy density. This analysis has been extended to SU(3) by C. Gong^^, who 
found the result 

hc^^g^iEp), (45) 

where (Ep) denotes the average energy per lattice plaquette. For a thermalized system 
the relation (Ep) = holds in SU(3); hence one obtains an entropy growth rate 
for SU(3) gauge theory of 

dS/dt = {h) ~ 0.54 g'^T. (46) 



This value coincides numerically with (twice) the thermal damping rate of long- 
wavelength gluons 

7i°^ = ^-^9'T:^ 0.264g^T. (47) 

The factor two between (46) and (47) may be understood by the remark that (47) 
gives the damping rate of the gluon amplitude; the decay rate of the probability is 
2jg. Although it is easy to see similarities between (h) and 2jg, their identity has 
not been established, and their numerical equality may be accidental. Independent 
of this, however, eq. (46) provides us with a precise value of the gluon thermalization 
time in the high-temperature, classical limit, defined as the characteristic entropy 
growth time near equilibrium: 



r« = {h)-' ~ 1.85/g^T. (48) 
Using the one-loop expression for the running thermal coupling constant 



one finds thermalization times of the order of 0.4 fm/c, as shown in Figure 17. This 

value is comfortably short on the time-scale of rclativistic heavy ion collisions, where 
the high-density phase is predicted to last for several fm/c. 
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Figure 17. Thermalization time of gluonic matter close to thermal equilibrium, as obtained from the 
Lyapunov exponent of thermal SU(3) lattice gauge theory*'^. 

The advantage of the definition (48) of the thermalization is that it is also mean- 
ingful for field configurations that are far from the average thermal configuration. 
Numerical studies of the rate of instability of coherent field configurations in SU(2) 
gauge theory have yielded even higher rates of entropy growth for the same energy 
density^''. More general investigations of this phenomenon would be helpful. 

FORMATION OF THE QUARK-GLUON PLASMA 

Approaches to the Thermalization Problem 

In the previous section we discussed how the quark-gluon plasma equilibrates 
when it is already close to thermal equilibrium. Relativistic heavy ion collisions pose 
a very different problem: how do the fully coherent parton wave functions of two 
nuclei in their ground states evolve into locally quasi-thermal distributions of partons 
as they are characteristic of the quark-gluon plasma state? There are mainly two 
approaches to this problem that have been extensively investigated: (a) QCD string 
breaking and (b) the partonic cascade. 



Figure 18. String-based picture of the formation of a quark-gluon plasma. Primary interactions lead 
to color flux tubes which break by quark pair production. 

In the string picture, developed from models of soft hadron-hadron interactions, 
one assumes that nuclei pass through each other at collider energies with only a small 
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rapidity loss (on average about one unit), drawing color flux tubes, or strings, between 
the "wounded" nucleons. If the area density of strings is low (not much greater than 
1 fm~^) they are supposed to fragment independently by quark pair production on 
a proper time scale of order 1 fm/c. Most realizations of this picture are based on 
the Lund string model^^, e.g. Fritiof^°, Attila^\ Spacer^^^ Venus'^^^ QGSM^"* and 
RQMD^^. When the density of strings grows further, at very high energy and for 
heavy nuclei, the formation of "color ropes" instead of elementary flux tubes has been 
postulated^^. 

Alternatively, a continuum description based on the Schwinger model of (1+1)- 
dimensional QED with heuristic back-reaction "chromohydrodynamics" -has been 
invoked to describe the formation of a locally equilibrated quark-gluon plasma'''''''*. 
One general aspect of these models is that initially part of the kinetic energy of the 
colliding nuclei is stored in coherent glue field configurations, which subsequently 
decay into quark pairs. The flux tubes carry no identifiable entropy. The entropy as- 
sociated with a thermal state is produced in the course of pair creation. In particular, 
there is no distinction between gluon and quark thermalization. 



Figure 19. Schematic view of a parton cascade (from ref. 85). 

The parton cascade approach, whose basic concepts were developed by BoaF^, 
Hwa and Kajantie^" and Blaizot and A. Mueller*^, is founded on the parton picture 
and renormalization-group improved perturbative QCD. Whereas the string picture 
runs into conceptual difiiculties at very high energy, when the string density becomes 
large, the parton cascade becomes invalid at lower energies, where most partonic 
scatterings occur at energies that are too low to be described by perturbative QCD. 
Nevertheless, the two descriptions may well be just two diff'erent formulations of the 
same physical processes, since there exists a remarkable similarity between the states 
produced by QCD bremsstrahlung and QCD flux tubes®^. 

Let us look at the "big picture", illustrated in Figure 20, where we distinguish 
three regimes in the evolution of an ultra-relativistic heavy-ion collision. Immediately 
after the Lorentz contracted nuclear "pancakes" have collided, scattered partons de- 
velop an incoherent identity and evolve into a quasithermal phase space distribution 
by free streaming separation of the longitudinal spectrum (a). Rescattering of these 
partons finally leads to the thermalization after a time of the order of 1 fm/c. The 
thermalized quark-gluon plasma then evolves according to the laws of relativistic hy- 
drodynamics (b), until it has cooled to the critical temperature Tc ~ 150-200 MeV, 
where it begins to hadronize (c). The physics governing the evolution during the 
three stages is very different. In this section we will concentrate mainly on the first, 
i.e. the preequilibrium state, because it is here where much recent progress has been 
made in our understanding. 
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Figure 20. One-dimensional space-time picture of the evolution of an ultrarelativistic nuclear colli- 
sion, distinguishing three regimes: (a) Pre-equilibrium state, (b) thermalized quark-gluon plasma, 
(c) hadronization phase, z denotes the beam axis. The dashed lines indicate the colliding nuclei. 

Parton Cascades 

As mentioned above, a fast parton of energy Ep penetrating a quark-gluon plasma 
loses about 1 GeV of energy in every 1 fm traversed. Its thermalization is therefore 
of the order 



which can be many tens or hundreds of fm/c for a past parton. How, then, is it 
possible that two nuclei colliding at 100 GcV/u (RHIC) or 3000 GcV/u (LHC) can 
deposit enough energy while penetrating each other to thermalize within 1 fm/c 
proper time? The resolution of this almost paradoxical situation has two aspects: 
(a) The emphasis is on proper time r. r = 1 fm/c corresponds to a very long time in 

the laboratory frame for matter produced at high rapidity. To wit, thermalization 

measured in the rest frame of the fast parton takes only 



where m* is the effective parton mass, and 3T is the average thermal energy of 
a parton. 

(b) Nucleons contain partons of all possible energies. Soft partons, i.e. partons 
with a; <C 1 mainly contribute to matter formed at central rapidity in the cm. 
frame, while hard partons, i.e. those with x ~ 0.3 — 1, are the source of matter 
at large rapidity in the cm. frame. It is quite illuminating to consider what 
would happen , if it were possible to instantaneously "shatter" the complete 
parton wave function of a fast nucleon. Indeed, one may ideally think of the 
collision of two nuclei at ultrarelativistic energies as an event during which only 
the phase coherence of the initial parton distributions is almost totally destroyed, 
but their longitudinal momentum distribution remains rather unchanged. Let us 
sec, therefore, what would happen if the nuclear collision acts as such a "phase 
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filter" 83. 



Intermezzo: Nuclear Collisions as "Phase Filter" 



We calculate the rapidity distributions of those quantities of the parton distribu- 
tion which are relevant to experimental data regarding the final state. One candidate 
is the entropy distribution, dS/dy, the other is the rapidity distribution of transverse 
energy, dET/dy. For the calculation of dS/dy and dEx/dy let us consider the phase 
space density of partons of flavor i with degeneracy rfj: 

Here A denotes the mass number of the colliding nuclei and V2 is the effective nuclear 
volume to be specified below. The momentum dependence of phase-space density can 
be derived from the measured parton structure functions: 

^ = i F,,A{xyr) « p ^^^^2^^ eM-pWo) (53) 

where we assumed the factorization of the longitudinal and transverse momenta: 
X is the Bjorken variable, the fraction of the nucleon's longitudinal momentum P 
carried by a parton. Fi{x) is the measured parton distribution in the nucleon, where 
shadowing effects are neglected here. One normally assumes a Gaussian distribution 
of width po ~ 0.3 GeV/c for the transverse parton momenta. 

To determine the spatial dependence of the phase space density, we use the 
concept of distributed contraction of partons^°. Let us call Y the beam rapidity 
and y the parton rapidity in the cm. system. Then only those partons occupy the 
Lorentz contracted region of fast moving nucleons with radius Rn which satisfy the 
uncertainty relation 

Here pl is the longitudinal momentum of the parton and Rn/ coshF = Rn ll is 
the longitudinal size of the contracted nucleons. Partons with smaller longitudinal 
momenta than (54) are more spread out in the longitudinal direction to satisfy the 
Heisenberg uncertainty relation, they form a cloud, containing mostly gluons and 
sea-quarks, surrounding the contracted nucleons (see Figure 21). 
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Figure 21. The longitudinal spread of the distributions of partons with rapidity y around a Lorentz 
contracted proton moving with rapidity Y . The spread is determined by the uncertainty relation. 

The effective spatial volume is then: 



(55) 



where y = miii(y,Y). One obtains a critical valnc Xc = cosh Y/{ARi^P) w 0.1. For 
X < Xc distributed contraction sets in. Considering the phase-space element 

W ^ W '^^^ ^' ^^^^ 
we can now calculate the rapidity distributions: 

VX pI cosh yY,di [fi In fi-{l± h) ln(l ± U)] , (57) 
VXPt cosh y^difi, (58) 
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where 

^Vpr^^^(-)^exp(-p|/2pg). (59) 

Using the approximation, SjFj(a;) D/x for sea partons, one obtains the fol- 
lowing expression for the integral (58): 

dEr ^ ADpo^^ 
dy tanhy 

The evaluation of dS/dy is more complicated and it is not expected to be linear 
in A, but it only depends on tanhy, as well. These results mean that the concept of 
''''phase filterincf^ leads to the prediction of a boost invariant phase-space distribution 
of partons in the rapidity region 1 < y < Y, where tanhy ~ 1, as it was postulated 
by Bjorken^'*. In the region < y < 1 we have a strong rapidity dependence in 
expressions (57,58), because the function (tanhy)"^ diverges in the limit y ^ 0. 
However, a detailed analysis^^ shows that this is balanced by the need to introduce 
a lower cut-ofF in the pT integrations, allowing for the continuation of the constant 
rapidity distributions down to y = 0. 
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Figure 22. Rapidity distributions of partons in collisions at RHIC (Au-|-Au at Scm (A=100 GeV), 
assuming complete "phase filtering" of tiie initial distributions, po parametrizes the initial transverse 
momentum distribution of partons in the nucleon. 



Let us consider a Au+Au collision at the characteristic energy of the Brookhaven 
RHIC machine (100 GeV/u), for which our considerations are likely to be applicable. 



As seen in Figure 22 an extended central plateau develops, where wc; predict a liadron 
multiplicity dN/dy > C~^{dS/dy) = 400, with C ~ 3.6 for a thermal ultra-relativistic 
Boss gas. This multiplicity may rise further due to rescattering and fragmentation. 
These effects can be studied systematically in a complete parton cascade model^^, as 
we will discuss further below. 

Parton Cascades (continued) 

Having discussed the initial parton structure of the colliding nuclei, let us exam- 
ine the scattering event more closely. When the two Lorentz-contracted nuclei collide, 
some of the partons will scatter and then continue to evolve incoherently from the 
remaining partons. Three aspects of these individual parton scatterings are worth 
discussing: 

(a) A parton-parton scattering can be described by perturbative QCD, if the mo- 

mcntiim transfer involved is sufficiently large. Nobody knows for sure where 
perturbative QCD becomes invalid, but typical choices*^"*^ for the momentum 
cut-off are p??'" ~ 1.7-2 GeV/c. 

(b) The elementary scattering has a finite space-time duration, which is given by the 
amount of off-shell propagation of exchanged virtual gluons or quarks. For the 
typical momentum transfer px = p™'" t^^is range is 

At, Ax ~ l/pf" ~ 0.1 fm(/c). (61) 

More precisely, this estimate applies to the QCD diagrams involving t-cliannel 
exchange, which dominate the total parton cross-section. For s-channel dia- 
grams, e.g. qq annihilation into a lepton pair, the elementary time scale is of 
order s^^/^, the inverse scattering energy in the parton cm. frame. 

(c) The time it takes for the scattered parton wave functions to decohere from the 
initial parton cloud depends on the transverse momenta of the scattered partons. 
Usually, one argues that the partons must have evolved at least one-half trans- 
verse wavelength At = t^/pt away from their original position before they can 
be considered as independent quanta. The considerations underlying this argu- 
ment are similar to those for the Landau-Pomeranchuk cffect^'^'^-'^. The critical 
issue here is that parton-parton cross sections are defined as squares of S-matrix 
elements which involve integration over all space and time from the infinite past 
to the infinite future. In the course of the nuclear interaction, however, that 
integration cannot be extended beyond their previous interaction point or their 
point of "formation" . If the plane wave factors representing the partons have not 
performed at least one-half complete oscillation the limited integral is not rea- 
sonably well represented by an S-matrix element. This is probably the weakest 
point in the line of arguments for a parton cascade and clearly would warrant 
more careful study. 

In any event, one concludes that a minimal time of the order of 0.1-0.3 fm/c must 
pass before scattered partons can be considered as incoherent field quanta, which 
fully contribute to the entropy and can rcscatter as independent particles. 

What happens after the initial scattering? The longitudinal parton momentum 
is little changed on average, i.e. ~ = xP, hence the scattered parton appears 
at rapidity 

y~ln(2p,/p^), (62) 

where p'rp is its final transverse momentum. Let us recall that the partons with 
momentum p^ were originally localized in a longitudinal interval (see Figure 20): 

1 1 



Hence, immediately after the nuc;lei have collided, the; fast partons are highly localized, 
whereas the distribution of soft partons is more fuzzy. Now let us make the outrageous 
assumption that the cloud of once scattered partons expands without any further 
interaction! Obviously, this implies the gradual separation of partons according to 
their velocity, because fast partons will leave the original interaction site quickly, while 
slower partons stay behind. Because of the approximate boost invariance already 
present in the initial state, we restrict the following discussion to the central rapidity 
interval, say, — 1 < y < 1. In particular all partons with momenta greater than \pz\ 
will have moved away from the central plane z = Q after a time 



Hence the longitudinal momentum spread of partons remaining at 2; = drops rapidly. 
It becomes equal to the average transverse momentum after: 

where we have assumed an average transverse momentum of 1 GeV/c, in line with 
results of numerical simulations ^^'^^ 

At that moment the distribution of scattered partons in the most central slab 
{z ~ 0) is approximately isotropic. It is also approximately thermal ? That depends 
on the actual phase space density of partons in the vicinity of 2; = 0. Only if every 
phase space cell up to \p\ ~ (p^) is occupied with probability close to one, can we 
speak of a quasi-thermal distribution. So let us estimate the phase space density. 
According to eq. (62), counting only the gluon distribution G{x), the number of 
partons in the original nuclear distributions, per unit of rapidity, is: 

^ ^ 2AxG{x) ^-=5= 6A. (66) 

ay d{m x) 

Here we have neglected nuclear shadowing, since we are aiming at a crude estimate. 

Let us further assume that the nuclei arc sufficiently "thick" and the collision energy 
is high enough that most available partons actually are scattered. The gluon phase 
space volume to be occupied is: 

= 167ri?l ^z — ,^ , ' , (67) 

where the factor 16 counts the color-spin degeneracy and Ra = tqA^^^ stands for the 
nuclear radius. Prom the uncertainty relation (63) we have Az ■ Ap^ 1; and since 
we have waited until Ap, has become equal to {p'-p), the rapidity spread of partons 
is Ay ~ 1. Combining everything, we find 

dVp/dy ~ -R\(p?). (68) 

TT 

Asking for full phase space occupation means requiring dVp/dy < dN/dy, or 



^1/6 > r_o{p^ ^ 2, 

TT 



(69) 



i.e. A > 50, if wc use (p^^,) ^ 1 GeV/c. Of eourse, not every available gluon from 
the initial structure functions will scatter when the nuclei collide. On the other 
hand, the fragmentation of scattered partons due to gluon bremsstrahlung effectively 
increases the number of scattered partons. Quantitative predictions can be made 
using a complete simulation of perturbative parton interactions, such as the parton 
cascade^^ or Hijing^'^. A calculation with Hijing for the system Au + Au at LHC 
energy was reported in ref. 88. The result was dN/dy ~ 800, {p'rp) ~ 1.75 GeV/c. 
The available phase space, according to eq. (68) is: 



dVp _ 2 /6fm- 1.75GeVY 
"rf^ ~ ^ 0.2GeVfm J 



1500, 



(70) 



i.e. phase space is about half occupied on average. The situation is less favorable with 
smaller nuclei or at lower energies. Figure 23 shows the number of primary parton- 
parton collisions for nucleon-nucleon collisions as well as for collisions of medium-sized 
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Figure 23. Area density of ijririiary jjartoii scatterings for various projectiles at RHIC and LHC 
energies, as a function of the reduced impact parameter b. Ra denotes the projectile radius. 



and heavy nuclei, respectively, in the RHIC and LHC energy ranges^^. The advantage 
of heavy nuclei, such as Au or Pb, for achieving a high density of scattered partons is 
obvious. Sulfur nuclei do not bring a great improvement over central nucleon-nucleon 
interactions. The increase in cm. energy between RHIC and LHC by a factor 30 
brings about a considerable increase in the parton density, as shown in Figure 24. 
The rise is not only due to the increased number of primary scatterings, but also 
partially caused by the growing amount of initial- and final-state gluon radiation. 
The calculations also show that only a small fraction of scattered partons are quarks. 
The assumption of fully thermalized quark distributions at very early times, r <C 1 
fm/c, is therefore most likely unwarranted*". 

Complete calculations following the evolution of the parton distributions micro- 
scopically until the attainment of thermal equilibrium have been carried out recently 
by K. Gciger*''''''-'^. A detailed account of these is given in Geiger's lecture at this 
School. Suffice it to mention here that he finds almost fully thermalized phase space 



distributions of gluons in An + Au collisions at RHIC energy (Ecm = 100 GeV/u) 
with T ~ 325 MeV after proper time r ~ 1.8 fm/c (see Figure 25). If this cap- 
tures the truth, there can be no doubt that a quark-gluon plasma will be observed in 
experiments at RHIC. 
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Figure 24. Distributions of scattered gluons predicted by HiJING for Au+Au collisions at RHIC 
and LHC. The predictions contains elementary QCD scattering cross sections and the radiative 
cascades of final state paxtons in the leading logarithmic approximation^^ . 



Figure 25. Rapidity and transverse momentum distribution of final-state partons as predicted by the 
parton cascade*^ ^ . The values of T axe obtained by a fit of the analytical isotropic fireball formula 
to the numerical results. 

Pre-Equilibrium Parton Physics 

The complete spectrum of phenomena occurring during the approach toward 

local thermal equilibrium still awaits exploration. Here I will discuss just two as- 
pects: color screening^^ and charm production^^, because they are dominated by the 



evolution of the giuon distribution, which is better understood at this moment. The 
experimentally very interesting contribution to lepton-pair production, which has also 
been studied^^, is more sensitive to the quark distribution. 

Color screenin^^: The screening length Aq of longitudinal color fields is important, 
because it defines the range over which coherent color-electric fields can extend. If 
it is larger than the characteristic confinement radius ~ 1 fm, long-range color- 
electric fields arrange themselves as flux tubes. In the one- loop approximation, Xd is 
determined by the momentum distribution /(k) of gluons according to: 

A-2 = hm / d'k^- q ■ Vk /(k). (71) 

A simple result can be obtained when the gluon distribution is exponential in kr 
and flat in rapidity y: 

2 _ 37rQs dN 



where dN/dy is the rapidity density of gluons. At the earliest time U ~ (fcr)"^, 
when the scattered gluons become incoherent, the calculations of ref. 88 based on 
the HiJiNG code predict 

AD^0.4fm (RHIC), 

Ar, w 0.15 fm (LHC). (73) 

We conclude that the color screening length at LHC energies will be so short, even 
after the first sequence of parton interactions, that coherent flux tubes cannot develop. 
The Lund model with its formation of independent strings that break up by creation 
of quark pairs is simply not applic;able under these conditions. At RHIC, on the 
other hand, the screening length appears to be marginally favoring a parton cascade 
description, but models with partially fusing flux tubes ("color ropes") may also be 
able to describe certain aspects of the pre-equilibrium phase of the nuclear collision. 

Charm production^^ : Charmed quark pairs are predominantly produced in collisions 
between two gluons. Since the charm production threshold is rather high, about 
3 GeV, the rate of production in a thermalized quark-gluon plasma with tempera- 
ture T ~ 300 MeV is negligible. Conventional wisdom until recently was, therefore, 
that most charmed quarks are produced in primary parton interactions^'*. In view 
of the high density of scattered partons with transverse momenta well above 1 GeV 
at collider energies, one may suspect that there is a sizable contribution to charm 
production from secondary parton collisions. This is borne out by a calculation of 
secondary charm production based on the initial scattered gluon distribution pre- 
dicted by HlJlNG^^. The additional charmed quarks populate predominantly the 
central rapidity plateau, where initial charm production is reduced by gluon shad- 
owing effects (see Figure 26). At LHC energies the total yield of secondary charmed 
quarks may be twice as large as that of primary charmed quarks. 

The amount of secondary charm production is sensitive to the thermalization 
time of the parton distribution. It depends on the ratio {cr c) / {'^tot) , where (ctot) is 
the average total parton-parton cross section that governs thermalization, while (ctc) 
denotes the averaged cross section for charm production, (crtot)"^ is proportional to 
the thermalization time rth- A measurement of the total yield of charmed particles 
(mostly D-mesons) in the central rapidity region would, therefore, provide valuable 
information on the time-scale of thermalization^^. 
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Figure 26. Rapidity distributions of charmed quarks produced in secondary gluon interactions, in 
comparison to the primary parton model predictions. The calculations are based on Duke-Owens 
structure functions with nuclear gluon shadowing. 



Plasma Evolution and Hadronization 

Once the; quark-ghion plasma has reached local thermal equilibrium, its further 
evolution can be described without reference to the parton reactions at the micro- 
scopic level. This concept was first quantitatively developed by Bjorken^^. The 
hydrodynamic equations for an ultrarelativistic plasma with P = admit a boost- 
invariant solution describing a longitudinally expanding fireball with constant rapidity 
density. When transverse expansion effects are taken into account, longitudinal boost 
invariance is partially destroyed, but the overall picture remains intact. This scenario 
has been thoroughly studied by a large number of theorists. Because the results are 
accessible in several fine reviews^^, I will refrain from discussing it here in detail. 

An important aspect of the late evolution of the quark-gluon plasma is its 
hadronization. Mostly it is assumed that the plasma expands and cools until it 
reaches the critical temperature Tc 200MeV and then converts into a hadronic gas 
while maintaining thermal and chemical equilibrium. More detailed descriptions of 
the dynamics of hadronization have been developed in connection with the problem 
of strangeness production, which is reviewed in the next chapter. 

A totally different approach, explored only very recently, consists in following 
the partonic reactions at a microscopic level, until the parton density has become 
sufficiently low to permit the formation of individual hadrons^^. A great deal is 
known about the mechanism of final state hadron production in e^e~- and NN- 
scattering but, unfortunately, we do not know whether this knowledge applies to the 
hadronization of a quark-gluon plasma in bulk. Therefore, the treatment of hadron 
formation at the end of a partonic cascade is frought with a great deal of uncertainty. 

At present, it is hard to judge the merits of either approach on a theoretical 
basis. Their validity and usefulness simply depends on whether the microscopic pro- 
cesses dtiring hadronization proceed approximately at thermodynamical equilibrium 
or not. The lesson from nuclear collision processes at much lower energies (below 1 



GcV/ii) has been that both scenarios are possible: ncnitron evaporation from highly 
excited nuclear fragments is well described by thermodynamics, but nuclear multi- 
fragmentation, where Coulomb forces play a dominant role, requires a more detailed 
theoretical treatment. 

QUARK-GLUON PLASMA SIGNATURES 

All theory of the quark-gluon plasma would be largely academic if there were no 
reliable signatures to observe its formation and to study its properties experimentally. 
It is impossible to present a complete review of quark-gluon plasms signatures here. 
I will, therefore, only try to capture the essential ideas and the current status of 
the theoretical studies on the most promising quark-gluon plasma signals. Anyone 
interested in more details is referred to the review of Kajantie and McLerran^^ and 
to the proceedings of the Strasbourg workshop^^. 

In order to shed some light on the connections between the many proposed quark- 
gluon plasma signatures I will group them in five categories, according to the physical 
properties of superdense hadronic matter to which they are sensitive. These are: 

1. thermodynamic variables measuring the equation of state; 

2. probes for chiral symmetry restoration; 

3. probes of the color response function; 

4. probes of the electromagnetic response function; 

5. "exotic" signatures of the quark-gluon plasma. 

Thermodynamic Variables 

The basic idea behind this class of signatures is to measure the equation of state 
of superdense hadronic matter, i.e. the dependence of energy density e, pressure 
P, and entropy density s on temperature T and baryochemical potential jiB- Here 
one wants to search for a rapid rise in the effective number of degrees of freedom, 
as expressed by the ratios e/T'^ or s/T^, over a small temperature region. These 
quantities would exhibit a discontinuity, if there were a first-order phase transition, 
and if we were dealing with systems of infinite extent. More realistically, we can 
expect a steep, step-like rise. According to recent lattice simulations this rise should 
occur over a temperature range of less than 10 MeV. 

Of course, one requires measurable observables that are related to the variables 
T, s, or e. It is customary to identify those with the average transverse momentum 
(pt), and with the rapidity distribution of hadron multiplicity dN/dy, or transverse 
energy dET/dy, respectively"*. One can then, in principle, invert the e — T diagram 
and plot {pt) as function of dN/dy or dEx/dy. If there occurs a rapid change in 
the effective number of degrees of freedom, one expects an S-shaped curve, as shown 
in Figure 27, whose essential characteristic feature is the saturation of (px) during 
the persistence of a mixed phase, later giving way to a second rise when the struc- 
tural change from color-singlet to colored constituents has been completed. Detailed 
numerical studies in the context of the hydrodynamical model have shown that this 
characteristic feature is rather weak in realistic models, unless rehadronization occurs 
like an explosive process^®. 

In order to trace this curve in nuclear collisions one probably has to vary the 
beam energy in rather small steps. This has not been done, so far, but it will be 
possible at RHIC. In nucleon-antinucleon collisions, however, one may make use of 
the existence of large fluctuations in the total multiplicity even from central N — N 
collisions. Using this tool, the E-735 collaboration at Fermilab^^" found a continued 



rise of {pt) for antiprotons and hyperons with multiplicity, reaching 1 GeV/c for 
the most violent events {dN/dy > 20) . When these data are analyzed in terms of a 
simple model, where one assumes that all hadrons are emitted from a longitudinally 
and transversely expanding fireball '^'^'^ , one finds that the surface velocity at high 
dN/dy must take on quite large values for the hadrons, reaching up to v/c = 0.8. 
Studying the hydrodynamical evolution that might lead to this final state, it is hard 
to believe that such a "flow" pattern can be produced at the level of hadrons, because 
the drag exerted by the dominant pions on the nucleons is far too weak to accelerate 
these to such speed. 
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Figure 27. Average transverse momentum of emitted hadrons as function of transverse energy 
d-Er /dy, representing the maximal energy density reached in a coDision. The different curves corre- 
spond to; (a) pion gas, (b) Hagedorn resonance gas, (c) quark-gluon plasma. 

The question then is: What produces the apparent transverse flow? Clearly, it 
must be established at the quark-parton level. It could be a consequence of expansion 
of a quark-gluon plasma or mixed phase. This could be tested by inspection of the 
full Pt spectra at high multiplicity^°^ . Alternatively, the transverse "flow" might be 
generated by the superposition of several extended minijets, as argued by Gyulassy 
and Wang^""*. It is not entirely clear that the two pictures are substantially different. 
Minijets might be the microscopic mechanism by which the transverse expansion of 
a quark-gluon plasma is produced. 

Models of the space-time dynamics of nuclear collisions need independent confir- 
mation, especially concerning the correctness of their geometrical assumptions. Such 
a check is provided by identical particle interferometry, e.g. of -k-k,KK, or NN 
correlations^'^'*, which yield information on the reaction geometry. By studying the 
two-particle correlation function in different directions of phase space, it is possible 
to obtain measurements of the transverse and longitudinal size, of the lifetime, and of 
flow patterns of the hadronic fireball at the moment where it breaks up into separate 
hadrons. The transverse sizes found in heavy ion collisions*"^, as well as in iV — A^ col- 
lisions at high multiplicity*"'^ are larger than the radius of the incident particle, clearly 
exposing the fact that produced hadrons rescatter before they are finally emitted. If 
interferometric size determinations would be possible on an event-by-event basis when 




Pb or Au beams become available, the correlation of global parameters like {pt) and 
dN/dy with the fireball geometry could be performed for each individual collision 
event. This will allow for much more precise study of the thermodynamic properties 
of superdense hadronic matter and may prove to be a sharp tool in the experimental 
search for a phase transition. 

Chiral Symmetry Restoration 

The two most often proposed signatures for a (partial) restoration of chiral sym- 
metry in dense hadronic matter are enhancements in strangeness and antibaryon 
production. The basic argument in both cases is the reduction in the threshold 
for production of strange hadrons (from about 700 MeV to 300 MeV) and baryon- 
antibaryon pairs (from about 2 GeV to almost zero). As Rafelski pointed out over a 
decade ago, the optimal signal is obtained by considering strange antibaryons, which 
combine both signatures^''®. The enhanced strange quark production in a chirally re- 
stored, deconfined quark-gluon plasma^*^^ leads to chemical equilibrium abundances 
for all strange hadrons, which would be difficult to understand on the basis of hadronic 
reactions alone^*^*. 

It has also been pointed out that strange particles, and especially antibaryons, 

would be produced more abundantly, if their masses would be modified even in the 
hadronic phase due to medium effects^°^. As discussed in an earlier section, the mass 
of K-mesons can be substantially lowered at finite baryon number density and the 
effective mass of antibaryons might be substantially reduced. 

An enhancement of strange particle production in nuclear collisions has been ob- 
served by many experiments^^°. However, we have also learned that such an enhance- 
ment alone does not make a reliable signature for the quark-gluon plasma. Strange 
particles, especially K and A can be copiously produced in hadronic reactions before 
the nuclear fireball reaches equilibrium. This mechanism seems to work very effi- 
ciently at AGS — as well as SPS — energies^^^'^^^. The cascade leading to enhanced K 
and A production has been studied in detail in the framework of the RQMD model, 
where it was found that most of the enhancement comes from ttN reactions initi- 
ated by the pions produced in first NN collisions, which have a strong nonthermal 
spectrum^ 

Maybe the most spectacular data in this respect are those obtained by the WA85 
collaboration^ at CERN, who find the following abundance ratios at mid-rapidity 
and for momenta Apr > 1 GeV/c: 



It is difficult to understand why the production of the doubly strange S should be 
particularly enhanced in a hadronic scenario, because the medium effects are expected 
to act mainly on the light quark content of baryons. Lattice calculations and sum 
rule estimates indicate that the light qq condensate is more rapidly depleted than the 
ss condensate in the medium. Rafelski^^** has argued that the observed ratios corre- 
spond to those found in a quark-gluon plasma that is about half equilibrated in its 
strangeness content. More recently, Rafelski and Tounsi"'^^''' have argued that strange 
baryon ratios seen by WA85 and other CERN heavy ion experiments can be consis- 
tently explained either by a quark-gluon plasma or hadronic gas with the parameters 
T = 220 MeV and = 340 MeV. Davidson et al."'^ have also pointed out that 
the ratios found by WA85 are close to those of a hadronic gas with effective volume 



A/A = 0.13 ±0.03 
S/A = 0.6±0.2, 



S/S = 0.39 ±0.07; 
S/A = 0.2 ±0.04. 



correction in complete chemical equilibrinm. However, even if this is so, the funda- 
mental question raised by the data is: how chemical equilibrium be attained during 
the short life of a hadronic fireball in any other way than through an intermediate 
quark-ghion plasma phase? 

Recently, attempts have been made to explain the enhanced A production seen 
by NA35 at midrapidity in terms of new mechanisms in the framework of collision 
models based on the string picture. Aichelin and Werner-'^^^ invoke the formation of 
"double strings" connected to the same leading quark to enhance the production of 
baryons containing strange quarks in the VENUS code. H. Sorge et al^^* introduced 
a mechanism for string fusion into "color ropes" which break faster and more often 
produce strange quarks and diquarks, into the RQMD model. This leads not only to 
strongly enhanced A production but also to a significant increase in the prediction 
for the number of produced antiprotons. Unfortunately no data on p production at 
the CERN-SPS are presently available to test this prediction. 

Strangeness enhancement has also been seen in the (^meson channel by the NA38 
experiment"'^'^^. Koch and Heinz^^" have argued that this effect can be understood 
as addition ^-production due to rescattering of secondaries, in combination with the 
small absorption cross section of the (/>-meson. The required density of scatterers 
agrees well with that invoked for explanation of the observed J/'^-suppression (see 
below). 

Color Response Function 

The basic aim in the detection of a color deconfinement phase transition is to 
measure changes in the color response function 



Although this correlator is not gauge invariant (except in the limit q ^ 0), its struc- 
ture can be probed in two ways (see Figure 28): 



Figure 28. Medium corrections to (a) the heavy quark-antiquark potential and (b) the energy loss 
of a penetrating parton (QCD jet) are sensitive to the color structure of dense hadronic matter. 

1. The screening length XdS'^'' = ng§(0)~^/^ leads to dissociation of bound states 
of a heavy quark pair, such as (cc). 

2. The energy loss dE/ dx of a quark jet in a dense medium is sensitive to an average 
of n°^(g^) over a wide range of q. 

Let me begin with the energy loss of a fast quark in the quark-gluon plasma, 
which was first studied by Bjorken^^^ in perturbative QCD. The connection between 
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energy loss of a quark and the eolor-dielectrie polarizability of the mecUimi was re- 
cently investigated by several authors^^'^'^ in analogy to the theory of electromagnetic 
energy loss. The basic formula is: 
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where jt denote the longitudinal and transverse components of the color-dielectric 
given explicitly in eris. (30). C is the Casimir operator for the color representation 
of the penetrating particle (C = | for quarks, C= 3 for gluons). Using eqs. (30) the 
expression can be evaluated analytically, yielding the energy loss of a fast quark; 
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This result includes both single-particle coUisional energy loss and energy loss through 
excitation of plasmons. All the medium dependence resides in the cut-off momenta 
q± in the logarithm: g+ w 2\/TE, while q- is a function of the Debye screening mass 
ruD = in the quark-gluon plasma phase. The magnitude of the energy loss is 
critically influenced by the strong coupling constant as, whose value unfortunately 
is not well known. The choices 0.2 and 0.3 preferred by various authors^'^^^^ leads to 
a stopping power between 0.4 and 1 GeV/fm for a fast quark. This may be a little 
smaller than the energy loss of a fast quark in nuclear matter. 

In addition, a fast quark loses energy by radiating gluons. Although this mecha- 
nism is strongly suppressed a high energy by the Landau-Pomeranchuk effect, it still 
contributes of the order of 1 GeV/£m to the energy loss^^. Adding the two contri- 
butions it thus appears that the stopping power of a fully established quark-gluon 
plasma is probably slightly higher than that of hadronic matter. However, in the 
vicinity of the deconfinement phase transition (if it exists!) there might be a region 
where the stopping power of strongly interacting matter decreases with growing en- 
ergy density. One would expect this effect to be particularly pronounced if the phase 
transition is of second order. The critical opalescence would in this case strongly 
suppress the emission of gluon radiation from the fast parton. If such an effect could 
be observed, e.g. by varying the transverse energy produced, it would clearly point 
toward a strong structural rearrangement in dense hadronic matter, possibly toward a 
quark-gluon plasma which is dominated by the propagation of collective color modes. 

The suppression of J/V' production, originally proposed by Matsui and Satz"'''*, is 
based on a simple, yet elegant idea: The ground state of (cc) pair does not exist when 
the color screening length \d = ^/gT is less than the boimd state radius ('"z/^)^^^ 
(see also Figure 12). Lattice simulations of SU(3) gauge theory^^"^'^^^ show that this 
condition should be satisfied slightly above the deconfinement temperature {T/Tc > 
1.2). The screening length appears to be even shorter, when dynamical fermions 
are included in the lattice simulations^^^ . In addition, the D-meson is expected to 
dissociate in the deconfined phase, lowering the energy threshold A.E* for thermal 
break-up of the J/i^. Blaschke^^^ has estimated, using the kinetic relation 

that the dissociation probability jumps significantly already at T^, and reaches unity 
at T/Tc ~ 1.2. 

The J /tp may still survive, if it escapes from the "dangerous" region before the 
cc pair has been spatially separated by more than the size of the bound state, i.e. 



more than about 0.5 fm. This may happcui, if the qiiark-giuon plasma cools very fast, 
or if the J/^f) has sufficiently high transverse momentum^^"^: pt > 3 GeV/c. The 
details of J/tp suppression near To are quite complicated and could require a rather 
long lifetime of the quark-gluon plasma state before becoming clearly visible^^*. 

On the other hand, the J/ tjj may also be destroyed in a hadronic scenario without 
phase transition by sufficiently energetic collisions with comoving hadrons^^®, leading 
to dissociation into a pair of D-mesons. This mechanism has recently been analyzed 
carefully by Gavin^'^" and by Vogt et al.^^^. In addition, the dependence of the 
suppression factor 5 on transverse momentum of the J/ tp is explained by a broadening 
of the transverse momentum distribution of pro jectile giuons due to prescattering"'^^^. 
Drell-Yan data^^"^ of Fermilab experiment E772 indicate that this proceeds like a 
random walk leading to a broadening A(p|,) which grows like These effects 

combine to explain most of the NA38 data. The result of these studies"'^^^'"'^'^'^ is that 
the pattern of J/tp suppression observed in experiment NA38 at CERN^^* can be 
understood on the basis of "standard" hadronic interactions, if one assumes comoving 
hadronic matter at density of at least 1/fm^ and an absorption cross section of the 
order of 2 mb. 

Karsch and Satz have recently analyzed whether this ambiguity persists at RHIC 
and LHC energies"'^^^. Assuming that a hadronic phase can be formed at energy 
densities of 4 — 7.5 GeV/fm'^, they find little difference in S{pt) in the accessible 
Pt range. However, they predict a substantial difference in the suppression of the T 
resonance at LHC energies^^^. 

Electromagnetic Response Function 

Electromagnetic signals for the quark-gluon plasma are in many respects ideal 
because they probe the earliest and hottest phase of the evolution of the fireball, 
and are not affected by final state interactions. Their drawbacks are (a) the rather 
small count rates and (b) the relatively large backgrounds from hadronic decay pro- 
cesses, especially tt" and r] decays. Electromagnetic signals probe the structure of the 
electromagnetic current response function: 



In the hadronic phase, n^i^(g^) is dominated by the p° resonance at 770 MeV, whereas 
perturbative QCD predicts a broad continuous spectrum above twice the thermal 
quark mass rriq = gT/^Q- At low <C 100 MeV collective modes are predicted 
to exist in both phases. In first approximation the collective quark-gluon plasma 
excitation, the plasmino^^"^ , has a somewhat higher effective mass than the collective 
7r"'"7r~ mode'^"^*, but its influence is hidden under strong nonresonant effects of soft 
QCD interactions in the plasma that cause a strong increase at low (see Figure 29). 
Unfortunately, these interesting modifications below the mass of the vector mesons 
will probably be overwhelmed by background from Dalitz pairs"'^''''. 

On the other hand, the production of lepton pairs with large invariant mass 
in the quark-gluon plasma phase may be sensitive to pre-equilibrium phenomena, 
e.g. collective plasma oscillations of large amplitude^*. Such oscillations are known 
to occur in the framework of the chromo-hydrodynamic model where the collision 
energy is first stored in a coherent color field which later breaks up into qq pairs. The 
ensuing collective flow could enhance the production of lepton pairs of high invariant 
mass. 
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Figure 29. Electromagnetic response function of dense hadronic matter (a) and quark-gluon plasma 
(b), as function of the invariant mass -y/g^ of a virtual photon. The hadronic response function is 
dominated by the neutral vector meson resonances p° and The dashed line in (b) shows the 
contribution from free quarks, while the solid line includes QCD interactions (as=0.3). Both cases 
exhibit collective modes at low g^, which are not shown here. 

The suggestion by Siemens and Chin^'*" that the disappearance of the p-mcson 
peak in the lepton pair mass spectrum would signal the deconfinement transition has 
recently been revived^*^. The basic idea is to utilize the fact that the qnark-gluon 
plasma phase should exhibit the higher temperature than the hadronic phase, and 
therefore lepton pairs from the quark-gluon plasma should dominate at high pT over 
those originating from hadronic processes. Unfortunately, the reasoning probably 
breaks down when one allows for collective transverse flow. Because of its larger 
mass, the p- meson spectrum is much more sensitive to the presence of flow than the 
quark spectrum in the quark-gluon plasma phase"^^^. 

Nonetheless, the lepton pairs from p-meson decay can be a very tiscful tool for 
probing the hadronic phase of the flreball. Heinz and Lee^^'^ have pointed out that 
the p-peak is expected to grow strongly relative to the to- and (jy peaks in the electron 
pair mass spectrum, if the fireball lives substantially longer than 2 fm/c. This occurs 
because of the short average lifetime of the p (1.3 fm/c), so that several generations 
of thermal p'^-mesons would contribute to the spectrum. In the limit of a very long- 
lived fireball the ratio of lepton pairs from p°- and uj- decays would approach the 
ratio of their leptonic decay widths (11:1). The p/w- ration can therefore serve as a 
fast "clock" for the fireball lifetime. 

The widths and positions of the p, oj, and peaks should also be sensitive 
to medium induced changes of the hadronic mass spectrum, especially to precur- 
sor phenomena associated with chiral symmetry restoration. This has been studied 
extensively The general conclusion, however, is that these modifications are prob- 
ably small except in the immediate vicinity of the phase transition. Changes are 
predicted to occur sooner, if the hadronic phase contains an appreciable net baryon 
density. E.g. a change in the K-meson mass could be detected via the induced change 
in the width of the ^-meson^^^. 

Direct photons, the second electromagnetic probe of dense matter, must face 
the formidable background from tt*^ and rj- decays. Whether these decays can be 
reconstructed and subtracted with sufficient reliability remains questionable, despite 
the remarkable achievements in this respect by experiment WA80 at CERN^^^. But 
even if it were possible, it is not clear what direct photons would tell. A new calcula- 
tion by photon emission from hadronic matter and quark-gluon plasma at the same 



temperature (T=200 MeV) by Kapusta, Lichard and Scibert^^^ has yielded virtually 
idential results for the two scenarios. This is by no means trivial, because the fireball 
is optically thin. The result is a consequence of the presence of thermal p-mesons and 
only occurs when the process up ti^ is taken into account. 

Altogether, the prospects for an unambiguous quark-gluon plasma signal from 
the electromagnetic sector are doubtful. Moreover, the count rates predicted at the 
future heavy ion colliders (RHIC and LHC) are quite small. The situation is reviewed 
in more detail by V. Ruuskanen^''^ . However, one should recall that electromagnetic 
probes are most sensitive to the earliest phase of a nuclear collision^°, and thus might 
prove to be a valuable tool for the detection of pre-equilibrium phenomena'^"' , even 
if they should eventually turn out not to be good probes of the thermal quark-gluon 
plasma. 

Exotica 

It would be nice if the formation of quark-gluon plasma would be associated with 

the appearance of completely novel phenomena: there would be no ambiguity in such 
signatures. Indeed one should remember that the proposal to look for quark-gluon 
plasma in nuclear collisions''*^ was originally derived from the apparent existence 
of unexplained phenomena observed in cosmic ray interaction, such as the famous 
"Centauro" events. The most probable exotic objects that might be formed from 
quark-gluon plasma are strangelets^^^^^ . As explained earlier, this name describes 
metastable objects with baryon number A > 2 that contain several strange quarks. 
The simplest such object is the strangeness S = —2 dibaryon, the if-particle, which is 
predicted to be metastable in the original MIT bag model'^' and might be produced 
in relativistic nuclear collisions'^^. Experiments'*^ searching for strangelets produced 
in relativistic heavy ion reactions are in progress at BNL, and in preparation at 
CERN. 

Recently there has been speculation about the possible formation of locally "dis- 
oriented" chiral vacua in relativistic nuclear collisions*'"''^. Such states would decay 
into a large number of pions, possibly with a strong isospin imbalance as observed in 
the Centauro events. They might be produced by some not well understood collec- 
tive emission process, or by spontaneous symmetry breaking when the dense hadronic 
system returns from the chirally restored phase. 

CONCLUSIONS 

It is appropriate to conclude this review of quark-gluon plasma physics by em- 
phasizing the positive aspects. Many of the proposed quark-gluon plasma signals 
have actually been observed already in the present experiments: J/ip suppression; 
enhanced production of strange hadrons, most notably of strange antibaryons; in- 
crease in transverse momenta of emitted particles. None of these results has been 
demonstrated to be an unambiguous signal of the quark-gluon plasma, so far. How- 
ever, one should bear in mind that the experiments were all performed with systems 
that were too small (Si and S are hardly "heavy" nuclei) and at energies too low 
to expect the formation of a full-fledged, sufficiently long-lived quark-gluon plasma 
state. In view of this, the experimental results are encouraging. 

On the theoretical side, a better understanding of what we really mean by a 
"quark-gluon plasma signature" is required'^^. In practical terms, we need a consis- 
tent formulation of what precisely is measured by J/^/i-suppressing and antihyperon 
enhancement. How do these signatures depend on the color response function or the 



quark c;orrdatioii function {qq) in the medium? Finally, we need to understand how 
pre-equilibrium processes influence predictions for the proposed signatures. These 
are difficult questions but, as I have tried to show, our understanding of the physics 
of the quark-gluon plasma has progressed to the point where these problems can be 
seriously addressed. 
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